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Abstract. A class of non-Hermitian quadratic su{2) Hamiltonians having an 
anti-linear symmetry is constructed. This is achieved by analysing the possible 
symmetries of such systems in terms of automorphisms of the algebra. In 
fact, different realisations for this type of symmetry are obtained, including 
the natural occurrence of charge conjugation together with parity and time 
reversal. Once specified the underlying anti-linear symmetry of the Hamiltonian, 
the former, if unbroken, leads to a purely real spectrum and the latter can be 
mapped to a Hermitian counterpart by, amongst other possibilities, a similarity 
transformation. Here, Lie-algebraic methods which were used to investigate 
the generalised Swanson Hamiltonian 1 are employed to identify the class 
of quadratic Hamiltonians that allow for such a mapping to the Hermitian 
counterpart. Whereas for the linear su{2) system every Hamiltonian of this 
type can be mapped to a Hermitian counterpart by a transformation which is 
itself an exponential of a linear combination of su{2) generators, the situation 
is more complicated for quadratic Hamiltonians. Therefore, the possibility of 
more elaborate similarity transformations, including quadratic exponents, is also 
explored in detail. The existence of finite dimensional representations for the su{2) 
Hamiltonian, as opposed to the su(l, 1) studied before, allows for comparison 
with explicit diagonalisation results for finite matrices. Finally, the similarity 
transformations constructed are compared with the analogue of Swanson's method 
for exact diagonalsation of the problem, establishing a simple relation between 
both approaches. 
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1. Introduction 

The interest in non-Hermitian Hamiltonians dates back to the early days of quantum 
mechanics. Traditionally these operators are associated with open systems, but, under 
additional conditions they can also be used to describe non-dissipative phenomena, 
e.g. Examples of such models have appeared in many areas, for instance in 

afhne Toda theories [4l[5] , Regge field theory [6] , atomic physics [7j , and also quantum 
spin chains |8j . In the mathematical community, non-Hermitian operators had already 
been studied [5HTT] , but the more recent work of Bender and Boettcher [T^] , connecting 
the reality of quantum spectra with an unbroken symmetry under combined parity 
and time reversal, led to a wide-spread interest in so-called T'T-symmetric quantum 
systems. Subsequently, it has been shown [T3] that unbroken PT-symmetry allows to 
introduce a new inner product with respect to which the time evolution generated by 
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the non-Hermitian T'T-symmetric Hamiltonian is unitary. The theoretical interest in 
non-Hermitian T'T-symmetric system also extends to classical physics, e.g. P^17) . 
Beyond mathematical physics motivations, recent activity in this area including 
experimental realisations of T^T-symmetric waveguides [18^ boosted the interest in 
the field. Possible applications are now envisaged ranging from optics jT9H2T] to 
condensed matter systems [22] and cold atoms f23 l l24 ) . 

The key feature of T^T-symmetric quantum mechanics is that real spectra and 
unitary evolution can be obtained for Hamiltonians not satisfying the Hermiticity 
condition. The connection between an unbroken T'T-symmctry and a real spectrum 
is a special case of a more general theorem for non-Hermitian Hamiltonians presenting 
any anti-linear symmetry |25II26) . If a Hamiltonian system possesses an anti- linear 
symmetry A, 

[i/,A]=0, (1) 
and its eigenvectors are also invariant under this symmetry, 

^ICn) = ICn), (2) 

we say that the symmetry is in its unbroken regime. Then, because (^|AU|^) = 

one can easily prove the reality of its spectrum without any reference to its Hermiticity 

properties [3S], P^ : 

^^nlCn) = H\Cn) = H A\C,) = AH\Cn) = A (S„|C«» = ^^^ICn)- (3) 

If a system satisfies ([T]) but not ^ , one says it is in a phase of broken symmetry. 
In this regime eigenvalues are not necessarily real anymore but, instead, come in 
complex conjugate pairs. U A widespread realisation in the community is to consider 
the anti-linear symmetry as the VT operator symbolizing simultaneous parity and 
time reversal transformations, together with complex conjugation. 

Nonetheless, having the reality of the spectrum established is not enough for a 
quantum theory. One still needs to guarantee its unitary evolution. If the system is 
non-Hermitian, 7^ H and a new Hermitian operator p, with respect to which the 
operator is Hermitian 0, 

H^p = pH, (4) 

must be included. This operator has the effect of changing the scalar product in 
the underlying vector space. Besides, it allows the construction of the 77 operator 
as p — rj, called the Dyson map, which makes possible the interpretation of the 
operators in the new Hilbert space, X and P, for instance, rather than x and p, as 
the genuine observables, 

H{x,p) ~ ri^^h{x,p) rj — h{rj^^x rj, rj~^prf) ~ h{X, P), (5) 

with the fundamental commutation relation [x,p] — [X,P] = ih preserved. In this 
case the non-Hermitian Hamiltonian H can be mapped to an isospectral Hermitian 
Hamiltonian via a similarity transformation 

h = rjHrj-'^ = h^. (6) 

According to the standard terminology, the Hamiltonian is said to be quasi- 
Hermitian if there is a Hermitian and positive metric operator p satisfying (|4]) , whereas 
for H to be considered pseudo- Hermitian p must be Hermitian and invertible [7ll27). 



f Note, however, that although the coalescence of the complex conjugate pairs into real values cannot 
be guaranteed, the latter might still occur, eventually. 
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The problem with the latter definition is that the definiteness of the metric is 
not guaranteed but with the former definition alone one cannot make a conclusive 
statement about the reality of the spectrum of the Hamiltonian. For more details 
about T^T-symmetry and pseudo- and quasi-Hermitian Hamiltonians see e.g. 
To circumvent these problems we impose simultaneously both requirements, leading 
to what will be henceforth referred to as pseudo-quasi-Hermiticity. 

Assuming that the r] operator has an exponential form, rj = for some operator 
T, the metric is automatically invertible and positive definite so that the reality of the 
spectra is guaranteed. The Dyson operator rj, therefore, constitutes a central element 
in the analysis of non-Hermitian systems. However, it is not unique and neither is the 
metric with respect to which the system is Hermitian according to (|4|). As a trivial 
example, one can multiply the metric by any unitary operator, as this leaves the 
Hermiticity of the Hamiltonian unaffected. Similarly, multiplying it by any operator 
that commutes with the Hamiltonian would yield a new metric. Nevertheless, we shall 
not be concerned with ambiguity issues for the moment. Instead, we will focus on a 
systematic method to construct the necessary elements for a fundamental description 
of the system, which is not trivial. There is in fact only a handful of known models for 
which the analytical non-approximative construction of an 77, or a metric, is possible. 
Often perturbative methods can be used to construct a metric for a non-Hermitian 
operator, e.g. [3T], but other approaches have been used, notably in [TJ[32H11|- 

The vast majority of investigations on PT-symmerty and non-Hermitian 
operators focusses on one-dimensional quantum systems whose Hamiltonians are 
composed of coordinate and momentum operators, that is, operators from a 
Heisenberg-Weyl algebra. There are some investigations on algebraic methods for 
sm(1, 1) |TJ|35] and recently also the E2 algebra ^43 . In Hermitian quantum 
mechanics an important class of Hamiltonians is composed of su{2) operators. These 
Hamiltonians have applications spanning diverse fields, reaching from the description 
of nuclear spins [44] to ultra-cold atoms (4^ . The spectral features of a PT-symmetric 
generalisation of a special example have been studied in [23) . 

In the present paper we investigate a general class of non-Hermitian Hamiltonians 
that are quadratic in the su{2) operators, and construct the corresponding 77 operators, 
in a continuation of work involving su{l, 1) quadratic Hamiltonians started in [T]. In 
the following section we report on the symmetries of the su{2) generators and present 
the possible anti-linear symmetries underlying su{2) Hamiltonians. As a by-product 
we can show in a unified manner that both su{l, 1) and su{2) algebras allow for three 
different symmetries, associated to the automorphisms of the algebra. We demonstrate 
that certain systems of interest present, for instance, a CPT-symmetry, yet to be 
defined, rather than just PT-symmetry. In the sequence, in section [3l we tackle 
the problem of finding rj for the quadratic Hamiltonian. We start with a simpler 
su(2) non-Hermitian Hamiltonian and construct metrics and isospectral Hermitian 
counterparts. Then, the analysis become more systematic as we investigate the use of 
metrics with exponents which are linear in the generators of the referred algebra, first 
to study linear Hamiltonians and later to consider the inclusion of quadratic terms. 
The following step is to explore more complicated similarity transformations capable of 
mapping non-Hermitian into Hermitian systems. In particular, we use exponents of the 
Dyson operator generated by elements of the enveloping algebra U {su{2)). Finally, the 
aforementioned results are compared with an explicit diagonalisation of Hamiltonians, 
thus allowing to make a correspondence between the pseudo-quasi- Hermitian approach 
and the Swanson method. 
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2. Symmetric Non-Hermitian Hamiltonians 

This section is devoted to exploring the symmetries ol Hamiltonians expressed in 
terms of certain operators. If the latter are generators of a particular Lie algebra, 
we investigate the transformations which leave the associated commutation relations 
invariant. For convenience, the algebras su(l,l) and su(2) are now studied, with 
respect to their automorphisms, in a unified way. 



2.1. A general framework for the invariant Hamiltonians 

In [T] an algebraic approach was implemented in order to find metrics for the 
quadratic generalisation of the so-called Swanson Hamiltonian. There, the method 
employed is based on the symmetry su(l, 1). However, one can easily adapt the same 
framework to su{2) algebra, whose generators can be realised in such a way that 
Bosc-Einstein condensates can be described [23]. Despite being very different from a 
representation theory point of view (the latter is compact, admiting finite dimensional 
representations, whereas the former noncompact), both algebras have a very similar 
structure, with three generators, Kq,K± and Lo,L±. In this section we will describe 
them in a unified way as much as possible and for that matter we introduce the 
M operators. Consider certain operators Mq, Mi, M2 satisfying the commutation 
relations 

[A/i, Afa] = i<jMo, [Mo, Mi] = iM2, [M2, Mq] = zMi. (7) 

Then, the operators Mo,M± — Mi ± iAl2 commute according to 

[A/o, M±] = ±Af±, [Af+, Af_] = 2crAfo. (8) 

It is clear that the choice ct = 1 corresponds to the su(2)-Lie algebra, with 
generators denoted by Mi = Li, whereas a = —1 represents su{l, 1), whose generators 
are taken to be Mi = Ki. There exist different automorphisms which leave these 
algebras, treated on the same footing, invariant |46| : 

Ti : Mo Mo, Af± -> M±, 

T2 : Ma -Mo, M± ^ -M^, (9) 

rg : Mo ^ -Mo, M± M^, 

and compositions of them. In fact, the mapping Hi — > —Hi and E±a — )• —E^ai of the 
Cartan subalgebra generators and step operators respectively, is an automorphism of 
any semi-simple Lie algebra due to the invariance of the root diagram under a — >■ — a. 
In this case it is an inner automorphism because the inversion is an element of the 
Weyl group. The first one is clearly just the identity transformation at the level of 
the M operators. However, it may still correspond to nontrivial symmetries, when 
expressed in terms of physical quantities in a particular representation of the algebra, 
as will be discussed below. 

The invariance of the algebras depicted above can be interpreted physically, in 
terms of coordinates and momenta, if we use representations in terms of the bosonic 
operators 



In cases where there is only one kind of boson, the index j may be dropped. Here the 
bosonic creation and annihilation operators , a can be defined as usual in terms of 
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coordinate and momentum coordinates x,p via (|10p. These operators, then, provide 
one with a representation of sm(1, 1) if one takes 

Ko = ^(^a^a+^y K+^^a^a\ K^^^aa. (11) 

In this representation we may reahze the automorphisms Q as 

at ^ -at, ^ -pr, 
at ia, <^ C, (12) 

at a, <=^ Txp. 

Thus, Tl is the usual VT transformation, an anti-hnear operation reversing 
time and space coordinates: x — > —a;, p ^ p, i — )• — z. Note that although this 
transformation can be achieved by the identity operator in (jH]), when writing the M 
operators in terms of bosonic operators, a nontrivial transformation associated to ri 
may be constructed with the use of a, at. The automorphism ri can be taken to be 
just the identity also in ((TU)) but the possibihty of using the equally valid VT operation 
is more interesting. 

A different transformation is introduced by T2, which changes the sign of the 
energy: uj -co. Having in mind the Dirac sea structure, one can interpret this 
operation as charge conjugation relating particles to antiparticles. One should not 
mistake this C with the one used in [55] , and denoted C , for a definition of the scalar 
product in the Hilbert space. The CPToperation can then be seen as composition 
of the first and second symmetries: CVT = tiT2 : Kq — —Ko,K± — —Kzp, a 
combination of charge conjugation with parity and time reversal transformations. 
The last transformation ra, on the other hand, is a symmetry in the phase space, 
intertwining coordinate and momentum: x — ^p, p — iujx. 

One may also represent the same algebra by means of to two distinguishable 
bosons 

Kq = ^ ^a\ai + ala2 + , = a\al, if_ = aia2. (13) 

The automorphism can then be described through 

t t 

Tl : ai — > — fli, a[ — > —a[, 



0,2 — > -02, 0,2 ^ ^02, 

r2 : ai — > a^, a\ — > 02, 

a2 — > —o\, aj — > — fli, 

T3 : ai ^ al, a| — > —02, 

a2 — > a\, aj — > — ai. 



(14) 



Again ti is the usual antilinear 7^7~-symmetry transformation, reversing 
coordinate and momenta of the same kind of bosons: i — > —i,Xi — > ~Xi,pi — > pi. 
In the coordinates and momenta representation, T2 can be achieved by considering 
LUi = u}2 = UJ together with i — —i,xi — ?> X2,Pi P2,X2 — >■ —Xi,p2 — ^ —pi- Also, the 
third automorphism, T3, relates position and momentum; Xi ±-^pi. 

Another similarity between the algebras su{l, 1) and su{2) is that the former can 
also be realised with the Schwinger Dyson transformation [47] 

Lq ~ -^{o\ai ~ ala2) , ~ a\a2, L_ = a^ai, (15) 

that maps the su{2) algebra to a two-dimensional Heisenberg-Weyl algebra with 
generators a\ 2 and 01,2- The most prominent example of this realisation are ensembles 
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of two-state boson systems, where creates a particle in state j, and aj annihilates one 
of them. Alternative realisations are given by two-dimensional single particle or one- 
dimensional two-particle quantum systems whose position and momentum operators 
Xj and pj are defined according to pUj) . 

In this representation the automorphisms satisfying ([9]) behave as 

Ti : fli — > — fli, a\ — > —a\, 

T2 ■■ oi ^ -ia2, al ial, , , 

. • t . • t y^^) 

02 — > Jfll, — > —iCl[, 

T3 : oi ^ 02, al al, 
02 — > ai, flj — 7- a|. 

In the ensemble of two-state systems realisation, the first symmetry is basically the 
identity, whereas the second and the third automorphisms correspond to interchanging 
the two states and can thus naturally be associated with parity. Associating the 
Heisenberg-Weyl algebra with position and momentum operators via (jlOp however, 
the automorphisms can be interpreted quite differently. 

Now, Ti still corresponds to a PT-symmetry (i — )■ —i,Xi — —Xi,pi — )■ pi) 
while the second automorphism T2 rotates the phase space by ^ in complex plane, 
wi -f-)- W2,xi -f-)- ±ix2,pi ±jp2- Finally, T3 just intertwines the particle's labels: 
wi O W2,xi O X2,xi ^ X2- This last one corresponds to the symmetry imposed 
in the system described in |23j . This realisation might be interpreted as involving 
bosons of different species transmuting into one another, bosons of the same kind 
but in different positions swapping their coordinates, or to describe a two-dimensional 
problem according to ([TU]). 

Although the algebras su{2) and su{l, 1) are not equivalent (since the former 
corresponds to a compact group and allows for a finite dimensional representation 
while the latter arises from a noncompact group and has only infinite dimensional 
representations), it is interesting to note that by applying Lq ^ — ^o, L± — > ^L^, 
which of course is not an automorphism, one can construct combinations which satisfy 
the commutation relation of the Ks, i.e. pT)) with a = —1. 

The framework described in this section is general and can be applied to other 
families of algebras in order to construct different anti-linear symmetries of the 
Hamiltonian. These symmetries, as we know, have an important role in selecting 
systems which potentially have real spectra and that consequently are interesting 
from a physical point of view. Moreover, it contains the whole class of Hamiltonians 
studied in this manuscript as a subcase and other families can be constructed similarly. 



2.2. Non-Hermitian su(2) Hamiltonians with an anti-linear symmetry 

Here we are interested in generic non-Hermitian Hamiltonians expressed in terms of 
su{2) generators, that is, Hamiltonians which are composed of angular momentum 
operators = Li, Ly = L2, = Lq, that fulfill the relations [L,j;,Ly] = iex,y,zLz 
with ex,y,z the Levi-Civita antisymmetric symbol. Alternatively one can use Lq = 
Lz, L± — Lx iLy, which commute as 

[io,i±]=±i±, [L+,L_]=2Lo. (17) 

The physical importance of su{2) Hamiltonians goes far beyond the description of 
actual angular momentum systems such as, e.g., nuclear spins [44]. The basis for 
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many other important physical reahsations of quantum su{2) systems is in terms 
of the bosonic operators pUj) . In Hermitian quantum mechanics quadratic su{2) 
Hamiltonians are of importance in fields such as nuclear physics and Bose-Einstein 
condensates [33] . These models are also very interesting due to the possibility of using 
semi-classical analysis in order to construct the scalar product of the non-Hermitian 
quantum system by relating the quantum Dyson operator with the Kahler metric for 
the corresponding classical Hamiltonian systems 24 . If possible, this could bring in an 
alternative approach to constructing metrics. Here we are interested in non-Hermitian 
generalisations of su{2) Hamiltonian systems. 

The most general complex Hamiltonian quadratic in the su(2) operators belongs 
to an 18 parameter family 

2 

H^^{aj+il33)Lj+ {(^jk+il3jk)LjLk. (18) 

j=a 0<j<k<2 

Any analytical function of Lj commutes with = + Ly + = Lq + ^{L^L_ + 
L_LjS). Thus the representation of such a function in the standard angular momentum 
basis has a block diagonal structure, where the n-th block is of size n x n. For the 
Hamiltonian (jl8|) . in particular, each block has a band structure. These Hamiltonians, 
however, are in general not pseudo-quasi-Hermitian. In what follows we will impose 
an additional anti-linear symmetry, which decreases the number of free parameters 
and allows for possible pseudo-quasi-Hermiticity. 

Below we comment on possible physical realisations and interpretations of the 
three automorphisms of the su{2) algebra. For the moment we shall investigate 
Hamiltonians fulfilling one particular symmetry as an example, t^, ©. Other cases 
can be treated analogously. Let us then consider Hamiltonians invariant under the 
third automorphism T3 combined with complex conjugation, that is 

Lq -Lq, L± izp, i -i, (19) 

or expressed in the angular momentum components 

Lx^Lx, Ly^Ly, Lz^-Lz, i-^-i. (20) 

Requiring the quadratic su{2) Hamiltonian (jl8p to be invariant under this 
transformation reduces the number of free parameters from 18 to 9. Explicitly, the 
Hamiltonians under investiagtion in the present paper are of the form 

H = iPoLa + (a+ + il3+)L+ + {a+ - il5+)L^ + 

+ aoo^o + (q^+o + il3+o)L+L(3 - (a+o - j^+o)i-io + (21) 
+ (a++ + i/3++)i+ + (a++ - iP++)L'i + a+^{L+L^ + L^L+), 
being also expressible in terms of ~ Lq, = ^{L+ + i_), Ly = 57 (i+ — L^): 

H = i'izLz + ^xLx + 'iyLy + 'izzL^z + ^xxL^ + '^yyLy + 

+ ^xy(LxLy -|- LyLx) + i^ixzLxLz + i^yzLyLz- (22) 
The exact matching between the two equivalent systems above is simply achieved with 
(ix = 2q!+, (,y = -2^+, qz = l3o, 

'ixx = 2(a+_ + a++), <iyy = 2(a+_ - a++), <;zz = aoo, (23) 
'ixy = —2/?++, i,xz — 2/3+0, '^yz ~ 2q;+o. 

In the linear part of the Hamiltonian (I21[) . the non-Hermiticity comes from the 
Lq term, differently from the su(l, 1) analogue. In that situation, the Kg term would 
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not break Hermiticity but the corresponding 'ladder' operators K± would, unlike here. 
Also, the fact that the Hermitian conjugate of L±Lq is not L^Lq, but rather LqL^, 
contributes for the breaking of Hermiticity in the quadratic part of (HI]). Clearly we 
can change the orders of these operators and gain extra terms in the direction of L± . 

Despite the fact of having H 7^ H\ we might still avoid complex eigenvalues. If 
the anti- linear symmetry of a Hamiltonian of this class is unbroken, it leads to a purely 
real spectrum and the existence of a Hermitian counterpart. For the construction 
of a subclass of non-Hermitian systems, the requirement of unbroken symmetry, 
however, is of little practical use, as it basically requires the explicit knowledge of 
the eigenfunctions. In what follows, we shall use a different approach to construct a 
subclass of non-Hermitian Hamiltonians, using Lie algebraic properties and an explicit 
ansatz for the 77 operator. This will allow us to construct not the whole, but a large, 
class of non-trivial pseudo-quasi- Hermitian Hamiltonians of the type (j2ip . 

3. Metrics and Hermitian counterparts 

Before we proceed to the general quadratic case, let us first illustrate the basic ideas for 
the analytically solvable case. In the study of an N particle two-mode Hubbard system, 
a natural Hamiltonian used to describe a Bose-Einstein condensate in a double-well 
potential emerges, which was made non-Hermitian in [2 3) by the introduction of an 
imaginary Lz term. It is instructive to represent the Hamiltonian in the components 
of the angular momentum operators as 

H = iq^Lz + ^xL^ + <;zzLl. (24) 

It was observed that in the vanishing interaction limit, with azz 0, the 
Hamiltonian can be expressed as 

H = _y^|T^e-arctanh(f-)L,^^garctanh(|-)L,^ ^2^^ 

indicating that it is in the same similarity class as the Hermitian Hamiltonian 
h = ^\/^x ~ ^z^zi with the operator rj — e^^^'^*^'^'^(~)^'', thus, acting as a Dyson 
map. The rj operator diverges at j^^l = |<^a;|. At this so-called exceptional point [481149] . 
the Hamiltonian is not diagonalisable, but has a Jordan block structure, where all 
eigenvalues and their corresponding eigenvectors degenerate. That is, at this point 
there is no equivalent Hermitian operator. For values of j^^j above this critical value, 77 
diagonalises the Hamiltonian again. For the interacting Hamiltonian, with a quadratic 
part, a metric operator of this simple form could not be found, though. 

While only some special systems can be solved via a direct diagonalisation, an 77 
operator of a simple exponential form that maps the Hamiltonian to a Hermitian, but 
not diagonal, counterpart can be found in certain circumstances. For this purpose, 
an explicit ansatz for rj is used to act adjointly on the Hamiltonian and the result 
is demanded to be Hermitian. This then yields constraints for the parameters of 
the Hamiltonian and of the 77 operator. Here we will investigate 77 operators whose 
exponent is linear in the s?i(2) operators and respect the anti-linear symmetry of the 
Hamiltonian. 

We start by analysing the linear Swanson-like model and then the more general 
families of quadratic Hamiltonians in terms of a metric with linear exponent with 
respect to the generators of the algebra. Later we introduce a metric with purely 
quadratic terms and finally we compare the similarity transformation underlying the 
pseudo-quasi-Hermitian approach with the Swanson method of diagonalization |33| 
and establish a connection between them. 
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3.1. Conditions for a linear exponent metric 

To find the whole class of quadratic Hamiltonians that can be mapped to a Hermitian 
counterpart by a simple t] with an exponent linear in Lj, an explicit ansatz for the 
Dyson operator is applied to the Hamiltonian and the result is demanded to be 
Hermitian. This then yields constraints for the parameters of the Hamiltonian and of 
the rj operator. In what follows we will investigate Dyson operators whose exponent 
is linear in the su{2) generators, that is, we make the ansatz 

?7 = exp[2(AoLo + A+L+ + A_L_)], (26) 

in terms the three parameters Ao,±. If we want the above transformation 77 not to 
break the eventual anti-linear symmetry of a Hermitian counterpart wave- function 1 0) 
for the non-Hermitian eigenstate |$) = 77^^ |(/)), we should impose 77 to be symmetric 
as well, implying Aq = iTo, A+ = A + iT, A_ = A — iV: 

77 = exp [2 {iToLo + (A + ir)L+ + (A - ir)L^)] 
^ exp{XxLx + XyLy + iXzLz), (27) 

with Xx — 4A, Aj^ — — 4r, A^ — 2Tq G R. The operator can be made Hermitian by 
choosing Fq = and its positivity is guaranteed by the condition 9"^ = 4(A^-|-r^)— r§ = 
j<j}^ = J {Xl + Xy — Xl) > 0. Due to the linearity of the exponent with respect to the 
su{2) generatores, it conserves the powers of the su{2) operators appearing in H. 

The adjoint action of the 77 operator on the L generators can be characterised 
by the matrix in the following equation: 

Lo 

L+ I , (28) 

whose determinant is simply unity and elements are given by 
600 = 1 + 8A+A_ , = - (2A^^)' , 6±± = (cosh0 ± Ao^)' , 

(29) 

6o± =T2A±5iifi« (cosh0±Ao^), b±o = T4.X^'-^ (coshO ± Xo'-^) . 

In some circumstances it is more convenient to use the metric ordered according 
to a Gauss decomposition such as 

^^^2k^L^^2k„Lo^2k+L+_ (30) 

In this case the parameters relate easily to the ones in (|26l) as \^ 

/ , , , smh0\ 2X+^ 2X-^ 

Kn — log COShy + An , K+ = . . „ , K_ — . , „ . (31) 

V " y' + cosh0 + Ao^ coshfl + Ao^ii^ 

We have already discussed that the non-Hermiticity of the Hamiltonian (PT|) 
results from the Lq term and from L^Lq and LqL_. However, when we impose 
h — h'' in we notice that that the adjoints of latter two terms can be recast 

§ This form is particularly interesting when acting on a highest weight state \fi) such that = 0. 

II Using the results in [l], the expressions above can be more easily calculated by noting that 
the they can be rewritten using a minor adjustment to reproduce now the algebra su(2) from 
1) generators according to Kq —Lq,K^ — > L^,K- —¥ — The comparison with the 
relations obtained previously for sm(1, 1) can be used to obtain the expressions above according to 
Ao — > ^ Aq, — ^ A — , A— — ^ ^'^+ Siiid kq — ^ — Ko,/^+ — >■ —kj^. 







( boo 


&0+ 


bo- 


L+ 


]-\ 


b+Q 


b++ 


b+ 


L- J 






b-+ 


6_ 



Metric operators for non-Hermitian quadratic su(2) Hamiltonians 



10 



in terms of the original operators after a simple commutation, introducing terms 
proportional to L+ and L_ . Therefore, all terms are Hermitian except for the linear 
components and there are three equations originated from these terms which need to 
be solved. After manipulation of the resulting constraints and the use of standard 
trigonometric identities, one can notice that they can be put in the form of fourth 
degree polynomial equations. For practical purposes, it is convenient to introduce 
the quantity Y = ISS^i^^ which is just an abbreviation for a natural combination 
appearing in the constraints. Because Y depends on in a transcendental way, solving 
for 9 explicitly remains a challenge. It is important to note that the quantity Y has no 
deeper physical meaning, with 9 (or ultimately, A,r,ro) being the important variable 
(or combination of variables) . These constraints are found to be of the form 

+ + + e^^V^ + y4 = (32) 

with i = 1,2,3 and coefficients given by the expressions bellow: 

= /3o, = 4r (2a+ - a+o) - 4A (2/3+ - /3+o) , 

^i'^ = 2To (2Aa+o + 2r/3+o - 4a+A - 4/3+r + /JoFq) , 
^i'^ = 4 (4 (r^ + A2) - r^) (r (a+o - 2a+) - A/3+0 + 2/3+A) , 

- (ro - 4 + A^)) (4ro (Aa^.o + r/3p,o - 2a+A - 2/3+r) + /3o (4 (F^ + A^) + T^)) , 

Co'^ = «+o, = -2 (ro/3+o + 2A (aoo - 2a+_ + 2a++) + 4r^++) , 

= 4 (6A (Aa+o + r/3+o) + To (F (aoo - 2 (a+- + 3a++)) + 6A/3++)) , 
= -4Aaoo (4 (r^ + A2) + T^) - 24ro (r^ + A^) /3+o - 2rlf3+o + 

+32 (r^A (a+_ + 3a++) + A^ (a+_ - a++) + r^(3++ - 3rX'^f3++) 

+8rl (Aa+_ + 3Aa++ + 3r/3++) , 

- 16 (r^ + A^) ((A - r)(r + A)a+o + 2rA/3+o) + 24rr2 (A/3+0 - ra+o) + (-a+o) 
+16ro (r^ (aoo - 2a+_ + 2a++) + TA^ (aoo - 2 (a+_ + 3a++)) + 

+2 (A^ - 3r2A) /3++) + 4Tl [T (aoo - 2a+_ + 2a++) - 2A/3++) , 

Co'^ = /3+0, = -4raoo + 2roa+o + 8r (a+_ + a++) - 8A/3++, 

= 24r (Aa+o + r^+o) - 4ro (A (aoo - 2a+_ + 6a++) + 6r/3++) , 
C^'^ = 2 (l2ro (r^ + a2) a+o + r^a+o - Sr^ (aoo - 2 (a+_ + a++)) 

— SrA^ (aoo — 2a+_ + 6a++) — 2rQr (aoo — 2a+_ + 6a++) 

+4A (3 (r2-4r2)+4A2) /?++), 

^ = 4 (SrA (r^ + A2) a+o + err^Aa+o - 4ro (r^A (aoo - 2a+_ + 6a++) 
+A3 (aoo - 2 (a+_ + a++)) + 2r^/3++ - 6rA2/3++) 

(-Aaoo + 2A (a+_ + a++) + 2r/3++)) - (24r2A2 + r;^ + 16A^ - 16r'') /3+o. 

Therefore, our strategy transformed the problem of determining the Dyson 
operator, the metric and the Hermitian counterpart associated to equation (pT|) . 
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ultimately, into that of solving the above set of coupled transcendental equations 
for the parameters a and /3 in the Hamiltonian, as well as A,r,ro introduced by 
the similarity transformation. Despite having completely different starting points and 
little in common, this technique can be seen as having the same spirit as the Bethe 
ansatz approach for integrable systems, for which the complete diagonalization of a 
problem is reduced to a set of coupled equations for the spectral parameters [50j . 

Once the Dyson operator is obtained the metric p — rf^rj can then be constructed 
as well as the Hermitian isospectral counterpart according to (jG]), providing us with 

-2A+^Lo + AooLl + 2A+^L+L^+ (33) 
+ {{A+ + Ao+ + i{B+ + Bo+)) L+ + {{A+ + Ao+ - i{B+ + Bo+)) i- + 
+ {A++ + iB++)LX + {A++ - iB++)Ll + 

+ (Ao+ + ^+0 + i{Ba+ + B+a)) L+Lo + {Ao+ + A+o - i{Bo+ + B+o)) LqL^, 

where each line in the right hand side of the equation above is Hermitian, given the 
coefficients Ai = cosh^ 9 Ai , Aij = cosh"' 9 Aij , Bj — cosh^ 9 Bj , Bij = cosh"* 9 Bij 
are real. In fact they can be expressed as 

A+ =a+{i- y2 (_4r2 + + ix^)) + 2y (/3o (r + ToXY) - p+ (ro + 4rAr)) , 

Am - "00 {Y^ (4 + A^) + r^) + l)' - BY (r/3+o + T^Y'' {AT^ (r^ (a+_ - a++) 
+ A2 (a+_ + Q'++) + 2TXP++) - A (4 (r^ + A^) + T^) /?+o) 
+ Y (-roA/3+o + 4 + A^) a+_ + AT^a++ - AX^a++ - 8rA/3++) 

+ y2 (r (4 (r2 + a^) + r^) /3+0 + sFq ((a - r)(r + a)/?++ - 2rAa++)) 

+ a+o (A + TToY) {Y^ (4 (r^ + A^) + Vl) + l)) , 

= 2Y {Y (rro (a+o (4y2 (r2 + A^) + l) - 16Aya++) 

+ TlY (A (a+o - 2Ayaoo + 4Aya++) - 2T'^Y [a^ + 2a++)) 

- (aoo - 2 (Aya+o + «++)) + rrgy^a+o - 8/3++ (A + rroy) (r - roAy)) 
+ /3+0 (r - ToAy) (y2 (4 (r^ + a^) + r^) + 1)) 

+ 2Ay (a+o + AX^Y^a+a - 2Ay (ago + 2a++)) 

+ a+_ (y4 (i6 (r2 + A^)' + r^) + 2tIy^ + 1) , 
A++ = a++ (y4 (16 (r^ - er^A^ + a*) + r^) - ergy^ + 1) 

- 2y (-r/?+o + y=^ (ro (a (-i2r2 + + 4a2) /3+0 

+ 2ro(A - r)(r + A) (aoo - 2a+_)) + 32rA (F^ - A^) /?++) 

+ y2 (r (-4r2 + Srg + 12A2) /?+o + SrroA (aoo - 2a+_) - 2TII3++) 

+ a+o (-i2rroA2y^ + a (1 - 3 {at^ + vl) y^) 
+ rroy ((4r2 + rg) y^ - 3) + ax^y"^) - 3roAy/3+o 

- 2y(A - r)(r + A) (aoo - 2a+_) + 2ro/3++) , 

io+ = 2y (2 (y (ro (A/3++ (4y2 (a^ - 3r2) + 3) - 4y (r^ + a^) /3+0) 

+ 4 (A^a+o + rA (/3+0 - 3Ay/3++) + T^Y 

+ rr^y (-4rya+o + 4Ay/3+o + + r^^A (-y^) /?++) 
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+ a+_ (A - TToY) (Y^ (4 (r^ + A^) + T^) + l) - r/3++ 

+ a++ (-i2rroA2y''' + a (3 {at^ + rl) - 1) 

+ aoo (rroF - A) (y2 (4 (r^ + a^) + r^) + 1)) , 
^+0 = 2Y (-ro/3+o + Y^ (4rA (4 (r^ + a2) + r^) /3+0 

+ 4ro (r^ (aoo - 2a+^ + 2a++) + TA^ (aoo - 2a+_ - 6a++) - 6r^Xp++ + 2X^/3++) 
+ r3(r(aoo-2a+_+2a++)-2A/3++)) 

- y2 (4r„ (r2 + A^) + rg/3+0 + 4 (r^A (aoo - 2 (a+_ + 3a++)) 
+ A^ (aoo - 2a+_ + 2a++) - 2r^f3++ + 6TX^13++) 

+ rl (A (aoo - 2 (a+_ + 3a++)) - 6T/3++)) + ATXYp+o 
+ ToY (r (aoo - 2a+- - 6a++) + 6A/3++) - 2Tf3++) 

+ a+o (Y"^ (leA^ - (4r2 + r^) 2) + SX^y^ + 1) - 2Ay (aoo - 2a+_ + 2a++) , 
B+ = 13+ (_4r2 -Tl + 4A2) + 1) - 2Ay {/3o + 4a+TY) + 2ToY (a+ + f3oTY) , 
B++ = 2Y (Fo (y (3r/3+o + Aa+o (-4y2 _ 3r2) - 3) + 4ry2 (r2 - 3A2) /?+o 

+ 4y(A-r)(r + A)(aoo-2a+_)) + 2a++) 

+ r (4Ay (aoo + y (8y(A - r)(r + A)a++ - 3r/3+o) - 2a+_) 

+ a+o (4y2 (r2 - 3A2) - 1)) + r2y2 (3A/3+0 + r (3a+o - 4Ay (aoo - 2a+_))) 

+ r3y2 (Aya+0 - ry/3+0 - 2a++)) + 8A3y3/3+o - 2XYP+0 

+ /?++ (y4 (r^ - 16 (r* - 6r2A2 + a^)) - 6r2y2 + 1) , 

Ba+ = -2Y (-2 (4y (F + ToAy) (a+o (A + TFoy) + /?+o (r - ToAy)) 

+ a+_ (r + ToAy) (y2 (4 (r2 + a2) + r2) + 1) 

+ a++ (-i2r2roAy^ + 4r3y2 + roAy (y2 (rg + 4a2) - 3) 

+ r (1 - 3y2 (r2 + 4a2)))) + aoo (r + roAy) (y2 (4 (r2 + a^) + r2) + 1) 

+ /?++ (-24rroA2y3 + a (6 (4r2 - rg) y2 + 2) 

+ 2rroy ((4r2-rg)y2 + 3) -8A'V2)) , 

B+o = 2y (rgy2 (a+o + y (-Aaoo + 2A (a+_ + a++) + 2r/3++)) 
+ r2y2 (6A/3++ - r (aoo - 4Aya+o - 2a+_ + 6a++)) 
+ 4ry (A (a+o + 4A2y2a+o - Ay (aoo - 2a+_ + 6a++)) 
+ r2y (2 (2Aya+o + a+_ + a++) - aoo)) + To (a+o {AY^ {T^ + A^) + l) 
+ 4y3 (-r2A (aoo - 2a+_ + 6a++) + A^ (2 (a+_ + a++) - aoo) - 2T^ (3++ + 6rA2;g+ 

- y (A (aoo - 2a+_ + 6a++) + 6r/3++)) + T (2 (a+_ + a++) - aoo) 

+ 2A/3++ (4y2 (a2 - sr^) - 1)) + /?+o (y* (- ((r^ + ax^) ^ - ler*)) + 8r2y2 + 1) , 

all subject to the constraining relations which will be determined in the following 
subsections. 

3.1.1. The subclass of linear Hamiltonians 



The system of nonlinear equations (j33p is not of trivial solution and we shall 
start by analysing the behaviour of the linear su{2) Hamiltonian according to this 
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framework. For that matter, we first consider the subclass of (132)) with a^^ ~ /3^^ = 
ct+o = P+o — 0, such that (pij) considerably simplifies to: 

H = il3oLo + (a+ + i/3+)i+ + (a+ - + aoo^o + i^+L^ + L^L+) (34) 

The non-Hermiticity of the operator above comes solely from the first term, all others, 
as a group, being invariant under Hermitian conjugation. A metric can be found by 
solving the constraining equations as long as «+_ = ^apo, meaning the above 
becomes just the linear Hamiltonian shifted by the Casimir operator L^. The Dyson 
operator, and consequently the metric, is then characterised by the parameters A, T, Fq 
such that the combination 9"^ — 4(A^ + F^) — Fq obeys 
tanh(6') _ /3o 



A - 4 [a+V - 13+X) 



(35) 



yi6(/3+A - a+TY + 8/3oFo(a+A + /3+F) - Pl{Tl + 4(A2 + F2)). 

The above constraints is to be used in (j33]) when calculating the Hermitian 
isospectral partner of together with the specifications of the linear Hamiltonian, 
namely a+-f = = a+o = /3+o = 0. A direct comparison between the non- 

Hermitian Hamiltonian and its Hermitian counterpart is presented in Figure 1, where 
we have specified all variables but one. The initial parameters in (p4|) are a+, a_ 
and /3o ■ By fixing the first two of them and reparameterizing the last one in terms of 
0, one can obtain a non-Hermitian Hamiltonian in terms of this latter variable. This 
choice is such that the free parameter measures the non-Hermiticity coming from the 
Lq term. Using a finite dimensional matrix representation of the generators of the 

type 

L-tjZ, to) = V "F 'rn){l im-f- l)|^,m± 1), LQ\l^m) — m\l,m) , (36) 

the diagonalisation of the Hamiltonian is possible and the plot of the eigenvalues 
as a function of 6*, is presented on the left. Here, it was used a representation of 
dimension 11 associated to I — h and m = — 5, — 4, 0, 4, 5. For the right hand 
side, the spectrum is obtained by diagonlizing the transformed Hamiltonian ([55)1 . 
with Q!+_|_ = /3+_|_ = a+o = /3+o = and subject to ([55)) . The diagonalization of 
this Hermitian counterpart involves other 3 free parameters, A,F,Fo, and 9 has an 
interpretation in terms of the parameters of the metric, so that the constraint ((35)) is 
always satisfied for 9"^ < 8. As one can see, the spectral curves for arbitrary values of 
the metric parameters, coincide, as expected, reinforcing the results presented. 

The subcase of ()34p where aoo = = corresponds to the non interactive 

limit of the Bose-Hubbard model when /3+ = 0. It is also the su{2) analogue of the 
so-called Swanson Hamiltonian [33 , defined for su(l, 1) 

Hs = fioKo + f^-K- + ^i+K+ (37) 

with energy spectrum — {n + ^) \j ^ — tJ'-t^+i in agreement with (fTT)) . Here ^o,± 
are nothing but parameters specifying an extension of the harmonic oscillator with 
characteristic frequency which can be recovered when = /i_ = 0; otherwise, 
when /i+ ^ non-Hermiticity arises. In a similar way as what was done for the 
Swanson model, one can express the Hamiltonian as an isospectral deformation of 
a multiple of the operator i^To, or Lq in this situation. The determination of this 
proportionality factor, similarly to equation (pS)). allows us to find the eigenvalues in 
the present case to be 

4+/3^)-f 



En=n\ {a\+Pl)-^, (38) 




Figure 1. On the left, the spectrum of the hnear non-Hermitian su{2) 
Hamiltonian I I34II with «_|_ = 4, /3+ = 4, qqo = Qh = held fixed and /3o = 

^'9^+16 ainh^^(fl)^^ ' "'^^'^'^ ^'^^ term responsible for the breaking of Hcrmiticity, 
being parameterized in terms of 8. An independent diagonalization is shown on 
the right, where we have the spectrum of II33I I. with 0++ = P++ = 0+0 = /3+0 = 
and subject to II35I I. again as a function of 9 in the interval [— 2v^, 2v^], with 

a+ = 4,/3+ = 4,000 = a+- = 0,/3o = gi + j^eslnh'^ie) besides A = 1, T = 
l,ro = V8 — 9'^. The graphs coincide, as expected, since the Hamiltonians are 
related by a similarity transformation characterized by A, F, Fq. 



in accordance with the result of the previous section and describing the behaviour 
of the spectra shown in Figure 1. The real energy spectra of these Hamiltonians 
are complementary: whereas g R for > 4^+^_ one has e M for 
jjQ = /3q < Ajl+jl^ = 4(0;^ + 0\) for the linear su{2) Hamiltonian. Alternatively to 
its direct derivation, the equation above could also be inferred by the correspondence 
between the generators L and K, mentioned in the end of subsection 12. II 

The Swanson linear Hamiltonian, however, is somewhat trivial and the quadratic 
terms are the ones giving rise to more interesting behaviour, such as interaction. 
Therefore the understanding of quadratic problems is relevant and shall be pursued 
in the next subsection. 

3.1.2. The general quadratic Hamiltonian 

Having presented the general framework allowing us to construct a consistent 
description of non-Hermitian Hamiltonians and shown the behaviour of the simpler 
linear subclass of Hamiltonians, we can now tackle the more involved situation where 
quadratic terms are added. To find the constraints for more general quadratic 
Hamiltonians pip to be mapped into a Hermitian counterpart by a Dyson operator 
of the form (|27p we again apply the ansatz to the Hamiltonian and demand the result 
to be Hermitian. 

It is important to note, that the rj operator preserves the powers of the angular 
momentum components appearing in the Hamiltonian. That is, it maps the linear 
part to a linear counterpart, and the quadratic part to a purely quadratic part, 
without linear terms. The quadratic part adds two additional constraints, besides 
those involving the coefficients of the terms linear Lq = in the Hamiltonian. 
Nevertheless, the fact that the Hermitian conjugate of LqL± is L^Lq implies that 
in order to rewrite the latter in terms of the former, the commutation relations need 
to used, introducing the non-Hermiticity of the quadratic parts in the linear ones. 
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As a consequence the constraint obtained in the linear case gets extended, with the 
quadratic part in general leading to additional condition, in the quadratic case but 
stays valid in the limit as a^j — >■ and Pij — >■ 0. Thus, although the introduction 
of the quadratic parameters increases the dimension of the manifold defined by the 
constraints, the range of linear coefficients for which a linear rj operator can be found, 
does not increase due to an additional quadratic terms, and lies within the region 

However, the class of ry operators (|27|) is larger than necessary. For any 
Hamiltonian that is mapped to a Hermitian counterpart by an 77 operator of the 
form (l39l) with Tq — 6 there exists another 77 with Fq = (5 = that also maps the 
non-Hermitian to a Hermitian Hamiltonian. This is because the Lq contribution for 
this particular kind of ansatz can be factored out, in a similar fashion as the Gauss 
decomposition of (j27p leading to pip , as a unitary operator, which cannot transform 
a non-Hermitian operator to a Hermitian one and vice-versa. Thus, to investigate the 
quadratic non-Hermitian Hamiltonian, we will use the simpler ansatz 

ij = exp [2 {(A -I- ir)L+ + (A - iF)i_}], . (39) 

where the Lq part has been removed. Further, it is useful to consider the parameters 
in the metric to be proportional, F — Xv. The exponent in the metric then reads 
T = 2A[(1 -|- + (1 — iiy)L^]. This rj operator can be interpreted as a rotation 

of the space x, y, z around an axis in the x, y plane by an imaginary angle, that is a 
boost. 

Here we present results for the Hermitian metric, with Fq = 0, but analogous 
expressions for Fq 7^ can clearly also be found, as argued before. We present results 
for the latter in an appendix section. By solving the constraints (|32p coming from the 
Hermiticity conditions on the transformed Hamiltonian ([6|), we end up with relations 
fixing the parameters in the metric, A and 

nu, through the quantity Q and also restraining some of the coefficients in the 
Hamiltonian. The metric transformations which map the non-Hermitian Hamiltonians 
into Hermitian counterparts must satisfy the following relation 

tanh(0) ^ (40) 
/3ovl + v"^ 



with = 2Av^^+T and further constraints 

{2a+ - a+o)(a+o - 2^+o('^ - 1)) 



2l3o{iy - 1) 



2/3+ - 13+0 + (a+o - 2a+)v "++ 1.2 _ ^ 

a+o(a+o - 4(/3+ - /?+o)) - (a+o + 2^+o(/3+o - 2^+))i^ - 2a+{a+o + 2/3+o - a+ov) 



2fio(v^ - 1) 

(/3+0 - a+av){{3+Q - 2/3+ -f (2q;+ - a.+a)v) - 4a++/3oi' 



2/3o(i^2 _ 1) 

Setting = in the expressions above, the constraints in the metric parameters 
and Hamiltonian coefficients simplify significantly, constituting simpler solutions. The 
linear limit (a^ = /3y = 0) reproduces as expected the results in (1351) . In these cases, 
for instance, the constraining equations might be more conveniently expressed in terms 



(41) 
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of tanh(2^?). If this is done, care must be taken in order to make sure no spurious 
solutions are added as there is an ambivalence between tanh(6') and coth(6') in this 
situation: tanh(26') = -r^^S£h|_ ^ ^cothe 

^ ' l+tanh^ 9 l+coth'^ 6 

The possibility of finding exact conditions for the so-called pseudo-quasi- 
Hermiticity to be present in a system is without a doubt a strong point of the 
current method. Moreover, this approach gives at the same time, not only information 
about the metric with respect to which the Hamiltonian has a unitary evolution with 
real eigenvalues, but also allows us to construct Hermitian isospectral partners. In 
order to do so, one just has to act with the Dyson similarity transformation on the 
original Hamiltonian, as described by dysonherm. The Hermitian counterpart for (I2ip 
constructed in this way is specified by p3p. with Fq = 0,r = vX, and constrained 
according to (|4T|). 

In Figure 2 below we confirm again that the original non-Hermitian quadratic 
su{2) Hamiltonian has a spectrum which exactly matches that of the constructed 
Hermitian counterpart with fj p9p . providing another consistency check. Similarly to 
the analysis of Figure 1, again we fix the parameters in order to satisfy the Hermiticity 
constraints (j40|) and leave one of them, /3o, free to vary. The diagonalization is done 
with the angular momentum representation for 1 = 5. However, now the intervals in 
both graphs are chosen to be different. Whereas the first graph gives a more general 
picture of the behaviour of the eigenvalues, the interval in the subsequent graph is 
also interesting as it shows in more detail that the apparent level-crossing of the 
characteristic values observed on the left actually corresponds to avoided crossing of 
levels, as shown on the right. Due to the presence of quadratic interacting terms, the 
structure of the energy levels is considerably richer than that of the linear case, in 
Figure 1, deviating from the simple harmonic oscillator. Differently from the linear 
case, now the Hamiltonian cannot be put in a form proportional to Lq only, so the 
analogue of equation ((38|) cannot be presented here. 

Whereas in the linear case the class of pseudo-quasi-Hermitian operators is 
identical to the class of operators for which an r] operator of the form ([27l) exists, and 
the latter can always be found via direct diagonalisation, these classes of Hamiltonians 
do not coincide in the quadratic case. Only a subclass of Hamiltonians quadratic in 
su{2) can be diagonalised by a similarity transformation with linear exponents in the 
generators of su(2). In addition there are non-Hermitian Hamiltonians of the type 
(|2ip for which such rj operators can be found, that map them to Hermitian, but not 
diagonal counterparts. A whole class of such Hamiltonians with a nontrivial spectrum 
can for example be obtained from the linear system, by adding a Hermitian quadratic 
part that commutes with the metric obtained for the linear part. 

It is also interesting to note that by fixing the parameters in the Hamiltonian 
arbitrarily and letting one of them, /3o again, vary freely, without respect to the 
constraints (|40|) we observe in Figure 3 the emergence of complex eigenvalues for the 
problem, a behaviour very similar to that described in |23| . This is a manifestation of 
the broken T'T-symmetry, in which case despite having a PT-symmetric Hamiltonian 
as we do, the eigenfunctions of the corresponding states at not invariant under parity 
and time reversal. For the values chosen a metric of the form ([26]) can only be found 
for either /3o = or /3o = 2.4733 but the graph shows us that the region of unbroken 
PT-symmetry is larger and we expect that a more general metric should exist to cover 
the whole sector. 

Considering the vanishing of the quadratic parameters a^, /3ij = 0, condtions 
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Figure 2. On the left, the spectrum of the hnear non- 
Hermitian su{2) Hamiltonian l|21|l with a-| = 1, «++ = 

1, a+o = 2, /3+0 = Ij Q!+ = 1 held fixed and /3+ = 
I ^2{2a+u - a+QU + /3+o) + Mu'^ + 1 (tanh (2Xyj'y'^ + l) + coth (2\y/v^ + l))] , 
which is the term responsible for the breaking of Hermiticity being parameterized 
in terms of 6. Independently, the diagonalization shown on the right depicts 
the spectrum of II33I I. again as a function of 9 in the interval [— 2\/2, 2\/2], 
with a+_ = 1, a++ = 1, a+o = 2, /3+o = Ij ol+ = 1 held fixed and /3+ = 
i \2{2a+u - a+oiy + /3+o) + /3o\/r^TT (tanh (2AVi^^ + l) + coth (2X^1^'^ + l))] 
together with the parameters of the metric r; to be A = 0.3, u = 0.2. Once 
more, the graphs coincide, a result that serves as a check that the similarity 
transformation constructed is correct. The different intervals used for the graphs 
shown give a general behaviour of the spectum, on the left, whereas on the right 
one has more details, particularly regarding the avoided level-crossing, which is 
unnoticed in the picture on the left. 




Figure 3. Imaginary part of the spectrum of quadratic non-Hermitian su{2) 

Hamiltonian II21I I as a function of /3o with a-\ = 1, Q++ = 1, Q+o = 2, /S+o = 

1, = 1, = 2. An equivalent spectrum for the Hermitian counterpart 1 133 II 
can be obtained, for instance, with A = 0.3, 1/ = 0.2 in fj. Due to the symmetry 
of the spectrum with respect to /3o, only a negative interval is shown. The energy 
levels cease to be real approximately when |/3o| > 3. 



IJ are automatically satisfied whereas can be rewritten as 



tanh (4VA2+T2 ) = 



(42) 



2(/3+A-a+r)' 

which coincides with psp . when Fq = 0, F = vX, written in terms of the double 
angle. Parametrising A = j(j)COS'd,r = ^(f>sind, the above equation can be 
interpreted as the orthogonality condition between the vectors (2a+, — 2/3+, /3o) and 
(sin 1} sinh </), cos sinh ip, cosh </)) [5T] . The second vector points from the origin to 
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a hyperboloid of two sheets, and can be written as {x,y,zty^l + x'^+ y"^) for any 
x^y G R. Thus, the absolute value of the slope of these vectors in the direction of 
z is always larger than 1. That is, if and only if the slope of (2a+, — 2/3_|_, /?o) in z- 
direction lies between —1 and +1 can it be orthogonal to a vector of this form. This 



gives the constraint |/3o| < 2^a\ + j3\, and it is now straightforwardly verified that 
the ry operator which was extracted from the diagonalisation fulfils this constraint, in 
agreement with p8|) . Due to the existence of more parameters in the general quadratic 
Hamiltonian, a simple geometrical interpretation of this form cannot be achieved. 

Now that a systematic study of the metric operators associated to the problem 
(pTt has been achieved, the constraints specified by ([iO]) and (|4T|) can be used to 
investigate the interacting Bose-Einstein Hamiltonian (|24|) , presented as a motivating 
example. This is simply accomplished by taking /3+ — = /3+o = ck+o — «++ = 

= in p3|) . As a result, the corresponding system ([33]) docs not have a consistent 

solution. This problem of not being able to find a metric of the proposed form 
motivates us to find alternative Dyson operators. 

Next, we shall discuss the existence of more general Dyson operators, by 
investigating the possibility of ry operators with an exponent including quadratic terms 
instead of only linear operators. 

3.2. Quadratic metrics 

There are examples found in [1] where one is faced with non-Hermitian Hamiltonians 
having real spectra but for which a metric of the form (|26p could not be found. This 
was only a confirmation of the fact that the ansatz used is far from being the most 
general. A priori there is no reason why a non-Hermitian quadratic su{2) Hamiltonian 
should be mapped to a Hermitian counterpart by an 77 operator of the simple linear 
form. From the known solvable examples of Heisenberg-Weyl type, it is for example 
well known, that the metric often is not an analytic function of the coordinate and 
momentum operators and might even involve nonlocal expressions. It is expected 
that the situation for su{2) Hamiltonians is in general not different. Thus, the class 
of non-Hermitian Hamiltonians of the type (|21|) with real eigenvalues is expected to 
be larger than the class fulfilling the constraints in the last section that is brought to 
a Hermitian counterpart by a simple 77 with a linear exponent. 

However, finding more general metrics soon goes beyond the capability of the 
approach followed up in the present paper, as we were using an explicit ansatz for the 
metric together with the Lie algebraic properties special to this ansatz. With that in 
mind we can choose an operator whose exponent is, instead, purely quadratic in the 
operators Lq, L±, 




gC-i-gCoioef+^+. 



(43) 



For that purpose we make use of the following actions 




(44) 
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The expressions above lead to 

rj L+ ri-^ ^ e'^° {cie'^'^"'"') , (45) 
rjL^rj-^ = e-^"e~'^^'"'"L- + 2(+e^" {cie^^""") + 

+ AC+e^<-'' (cie^^"'^'') Co - 4C^e'^° {cie^^°'"'f , 

where we have abbreviated 

co = Lo + 2C-L2_, (46) 
ci = L+- 2C-L- - AC-UL- - ACiLt. 

From this actions we can aheady extract a number of interesting observations. 
First of all, we observe that the first equation in (|331), when read from right to left 
so to speak, is already an example of a non-Hermitian Hamiltonian that is mapped 
to a Hermitian counterpart by a metric with a quadratic exponent. This is a trivial 
example, because in the standard basis of angular momentum it is immediately obvious 
that the eigenvalues of the non-Hermitian operator Lo + 2(^±L\ are identical to those 
of Lq, and thus real, as in this basis, it corresponds to an upper/lower triangular 
matrix with the eigenvalues of Lq on the diagonal. This example, is not of the type 
(|2ip we have been investigating in the present paper, since there we cannot vanish L_|_ 
or L_ independently. However, we can easily construct an example of this type, by a 
rotation, so that the Hamiltonian on the right side of the first equation in (|^^ is in 
fact of the type (PT|l . 

Secondly, due to the fact that an exponential of a product of su{2) operators does 
not preserve the power of the generators L upon which they act adjointly, the mapping 
above is clearly not a linear transformation. This strongly complicates the structure 
of constraints obtained from demanding the transformed Hamiltonian be Hermitian. 
In order to overcome this difficulty in formulating condition ([5]), another possibility 
is to consider for simplicity a matrix representation for the problem in an attempt to 
construct a metric with quadratic exponents, as in (j43l) . When calculating it, one can 
make use of the fact that, unlike su{\, 1), the algebra su{2) admits finite dimensional 
representation in order to analyse the system more easily by working with a specific 
realisation. The simplest of these is the two-dimensional spin ^ representation, 
generated by Pauli cr-matrices. However, the simplicity of such representation does 
not always play in our favour. If we make use of well known properties 

CTj • a J = (5ijl + i ^ e.yfeO-fc , (47) 

k 

we can immediately see that the inclusion of quadratic terms LiLj of the enveloping 
algebra U{su{2,)) in the Hamiltonian does not bring anything new to the problem 
when formulated in terms of 2 x 2 matrices. This is because the generators of su{2) 
correspond to a basis in two-dimensions but this property does not hold anymore 
for higher dimensions (higher spin representation). It seems that the higher the spin 
representation the more effect the quadratic terms will have. In the basis of the spin 1 
representation, with / = l,m = — 1,0, lin the expressions in take the simple 
form of 3 X 3 matrices, so that one is left with 9 equations, complex in general, when 
imposing the condition (|6]) . Requiring the quadratic metric (|43p to be invariant under 
the same anti- linear symmetry we must have C± = Ci ± iQ2 , with Co i Ci i C2 € K-? which 
makes the metric also Hermitian. 
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Hermiticity can then be imposed by requiring 



/3o 



4/3++C+ (8Cl - ^2 - 4C^ + 2) 



16C+ + 1 
/3+0 (2C+ - 1) 



-2C+ + e2Co(c_(4C+-2) + l)-l' 
a+Q (2C+ + 1) 



2^ 



(C- (Mi 



C- 



V2C+^/l-2C2 -C+ 



4Ci-i 



(48) 



with IC+I < Ahhough the actions (|^5|) used are representation independent, the 



Hermiticity constraints above (|48l) were constructed from a 3 x 3 representation and 
therefore are not vahd in generaL It is interesting to note that no situation for which 
the quadratic metric coexists with the hnear metric could be found. This reinforces 
the need of more general metric operators as the limitations of the Dyson operator 
with only linear operators in the exponent becomes clear now. Finally, in Figure 4 
we have an example of non-Hermitian Hamiltonian with real eigenvalues for which a 
quadratic metric, of the type (j43p . can be found. 




Figure 4. First levels of the real energy spectrum, as a function of /3o, of 
the non-Hermitian Hamiltonian H = i/SoLo + a(i + ^^^) {L^ - L^) + -^'o + 
3(-L-L+ + L+L—). This system can be mapped into a Hermitian counterpart 

with a Dyson operator of the form q = e e -e'.o^oe-i + and has eigenvalues 

j81/3g+64(431-160V7) 
12 and 7 it — — — gy= , therefore always real. 



3.3. Swanson diagonalisation 

It is in general not clear why a quadratic non-Hermitian Hamiltonian should allow for 
a simple 77 operator of the form (|26p at all, and we now know that there is a wider 
class of Hamiltonians for which only a more complicated 77 operator can be found. In 
what follows we shall briefly outline how the class of diagonalizable operators can be 
constructed, and then we will focus on the more general class of operators that have 
an r/ operator of the linear form, that does not necessarily diagonalise them. 
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For the investigation of the quadratic Hamiltonians that can be diagonahsed by 
the rj operator with hnear exponent, it is convenient to view the action of as a 
generahsed Bogohubov transformation If we make in (pS)) the identification 



= 2A_ 



(5 = 2A- 



sinh 9 



sinh 6 



13 cosh 9 + Xa 
7 = cosh 6* — Aq 



sinh 9 



sinh( 



(49) 



it is also possible to interpret the same transformation as a generalized Bogoliubov 
transformation 1331: 





\ j jf3 + da 


—aj3 7(5 


L+ 


^26/3 




L- / 


' \ 2ja 


-a" 7^ 



L_ 



(50) 



with — a5 — 1 io guarantee unity determinant, reducing the number of parameters 
in the transformation to 3, the same as in the similarity transformation. 

For the su{\, 1) representation in terms of Ki the interesting feature which allows 
the diagonalisation of the system |1I33| is that new creation and annihilation operators. 
For the representation in terms of L the corresponding new bosons are 



(51) 



As a result, the transformed operators are written in a very simple form in terms 
of these new bosonic operators 







( 7 


—a 







1 


ai \ 






-d 


/3 












Cl 










/? 


b 




A 


V C2 j 




I 





a 


7 } 


V 


4 / 



02^2), 



C2d\ 



(52) 



with [c?i, Cl] = 1, [(i2, C2] ~ 1 and the remaining commutators vanishing. Consequently, 
TSSi = Cidi behave as number operators and [-/Vi,Ci] = Cj, [A'^^, d^] — — ( 
Hamiltonian can be expressed as 



a ^L^La 



If a 



(53) 



?(-r 



H — OoLq + UooLq 

N2) + ^{Ni - 2N1N2 + Nl) + 

+ d+-{Ni + N1N2) + a-+{N2 + N1N2), 

then, despite being still non-Hcrmitian (iV/ 7^ Ni), it is already diagonahsed for some 
states |ni) €5 \h2) satisfying Ni\ni) ~ ni\ni). Equivalently ^23 the states can be 
represented in the usual angular momentum basis \l,m) with I = ni + n2 and m = 

n-i — n2, such that Lo\l, m) — m\l, rh) and L±\l, m) = \J (f=p rri)(l ± fin + l)|f, rh ± 1). 

The conditions for which the transformation brings the Hamiltonian into this 
diagonal form can be obtained from [1] having in mind that Kq — > —Lq, — > 
L^,K^ — > — i+. For the linear case, for example, the constraint obtained is almost 
the same as for the similarity transformation (1351) : 

tanh(6') _ a+\ + (3+T 
' 9 ^ a+V - P+X 
so that when Fq = one must have a+A + /3+F = 0, a condition which can be used 
to fix the ambiguities in the metric. In this case, the metric will also diagonalise the 
Hamiltonian. 



F, 
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4. Conclusion and Outlook 

In this paper we carried out the construction of metrics and Hermitian counterparts 
associated to a certain class of quadratic su{2) non-Hermitian Hamiltonians by using 
a Lie algebraic formulation. A method first employed for systems with infinite 
dimensional representation is now applied to problems admitting finite dimensional 
representations. In this way our methods and results can be successfully checked with 
explicit diagonalisation. The symmetry present in such systems is not VT but still 
it was possible to construct non-Hermitian models with real spectra by determining 
the metric operator and isospectral Hermitian counterparts. In fact, examples of 
non-T^T-symmetric models with real eingenvalues can be found in e.g. |52H54| but 
the construction of Dyson operators, as done in this manuscript, corresponds to an 
important step towards a complete understanding of the non-Hermitian structure of 
the Hamiltonian. 

The equivalence of the spectra of the two Hamiltonian partners, the original 
non-Hermitian and the transformed Hermitian one, was shown in comparative graphs 
encoding the algebraic structure determined in exact form. We have taken our metric 
ansatz beyond the linear operators in the exponent and explored the possibility of 
having a more complicated ansatz, quadratic in the su{2) generators. The difficulties 
imposed by this metric regarding the commutation relations of quadratic terms were 
overcome with the use of a finite-dimensional representation. We have used studied a 
3x3 example which has the advantage of being simple enough to allow for a solution 
but at the cost of not necessarily having a generic representation independent result. 

Finally, we have emphasised that the anti-linear symmetry of the Hamiltonians 
investigated in this manuscript correspond just to one of three possibilities presented 
and there are other different non-Hermitian Hamiltonians expressed as linear and 
quadratic combinations of su{2) generators which can have real spectra. The other 
possible realizations of such a symmetry still remains to be investigated in more detail, 
in particular the appealing situation where T'T-symmetry comes naturally combined 
with charge conjugation. Metrics and Hermitian counterparts for non-Hermitian 
Hamiltonians described by both su{2) and su(l,l) generators may also be treated 
with the knowledge of the transformations presented here and in fT , in combinations 
which are extensions of Jaynes-Cummings models |55[l56j . 

Appendix A. Metrics with Fq 7^ 
Appendix A.l. F = 

Letting aside the condition that 77^ = rj, we can still impose some simplifications by 
vanishing one of the parameters in the metric. Again we introduce some proportional 
factor relating the remaining parameters, Fq = Xv, such that T = 2X[ih'Lo + L^ + L^]. 
The solutions in this case are expressed as 



Y 



tanh(e') _ 4(/3+o - 2/3+) ± A 



(A.l) 



2{{2a+ ~ a+o)!^ - Poiiy^ + 4}) 



A 



Vl6(/3+o - 2/?+)2 + 4/3o(2a+ - a+o)u - l3o{u^ + 4) 



with 9 = AVI 



— v 



'2 and 



2a+- + — — + 



+ 



8Y 2 
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(A.2) 



P++ 8(i-r2i,2) + 8r(i-y2i/2) 

There is also another set of solutions for which one cannot take the y limit 
tanh(6») = ±i\/4-z/2 



_ ^ (a+o±^+o)(2+.-) (A-4) 
a _ _|_ ('aoo ^ ^ (a+oT;3+o)(4+12i/^+i/-') 

Appendix A.2. A = 

For this other two-parameters metric which remains to be analyzed we use Fq = Tv, 
and the argument in the metric becomes T = 2iF[z/Lo + — -^-]• 
tanh(£ ^ ^^(a+o^2a+)^ 
~ 2(4(/3+o-2/3+)z. + /3o(t^2+4)) ^-^^ 

A = Vl6(a+o - 2a+)2 - 4/3o(4(/3+o - 2/3+)i/ + Poi^^^ + 4)) 



with 6* = rV4 - j/2 and 

^ ^ /3+0 , g+o^ , (l + 2y2)(ct^p_;3 )_8ya++ 
aoo =2(a+_-a++)- — + ^ + 4y(y^ - 1) + 

(l + 2y2)(a+o + /3+o)+8ya++ 4Y{P+o +Ya+oiy) 

4y(yi. + i) i + y2(4 + z/2) ■ > 

_ iy{16Ya++ + /3+o(4 + ^2)) a+o(l + y2(4 + u^)) 

'^++ 8(y2!y2_i) 8y(y2z/2-i) ^ 

There is also another set of solutions for which one cannot take the u ^ limit 
tanh(6') = ±iV4- 1^2 



a+ = 
a++ 



— ±yPo-— 2D — ) 



(a+oT/3+o)(2+.^") (A-^ 
« _ ^ /aoo _ ^ _ (a+o±;3+o)(4+12i/"+i/*) \ 
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